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ABSTRACT 

Context. The precision of radial velocity (RV) measurements to detect indirectly planetary companions of nearby stars has improved 
to enable the discovery of extrasolar planets in the Neptune and Super-Earth mass range. Detections of extremely low mass planets, 
even as small as 1 Earth mass or below, in short-period orbits now appears conceivable in ongoing RV planet searches. Discoveries of 
these Earth-like planets by means of ground-based RV programs will help to determine the parameter //g,, the frequency of potentially 
habitable planets around other stars. 

Aims. In search of low-mass planetary companions we monitored Proxima Centauri (M5V) as part of our M dwarf program. In the 
absence of a significant detection, we use these data to demonstrate the general capability of the RV method in finding terrestrial 
planets. For late M dwarfs the classic liquid surface water habitable zone (HZ) is located close to the star, in which circumstances the 
RV method is most effective. We want to demonstrate that late M dwarfs are ideal targets for the search of terrestrial planets with the 
RV technique. 

Methods. Using the iodine cell technique we obtained differential RV measurements of Proxima Cen over a time span of 7 years 
with the UVES spectrograph at the ESO VLT. We determine upper limits to the masses of companions in circular orbits by means of 
numerical simulations. 

Results. TheRV data of Proxima Cen have a totalrms scatter of 3.1 ms"' and a period search does not reveal any significant signals. 
In contrast to our earlier results for Barnard's star, the RV results for the active M dwarf Proxima Cen are only weakly correlated with 
Ha line index measurements. As a result of our companion limit calculations, we find that we successfully recover all test signals with 
RV amplitudes corresponding to planets with m sin i > 2 - 3 Me residing inside the HZ of Proxima Cen with a statistical significance 
of > 99%. Over the same period range, we can recover 50% of the test planets with masses of m sin i > 1 .5 - 2.5 M^. Based on our 
simulations, we exclude the presence of any planet in a circular orbit with m sin / > 1 MNcptune at separations of a < 1 AU. 

Key words, stars: individual: Proxima Cen - stars: planetary systems - techniques: radial velocities 
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envisioned for the formation of terrestrial planets around main I" this paper we present 7 years of high precision RV data for 

our closest neighbor in space: the M5V star Proxima Centauri. 

Send ofprint requests to: M. Endl We demonstrate that with the data in hand we could have already 

* Based on observations collected at the European Southern detected planets with minimum masses as small as 1 - 2 M®. 

Observatory, Pai'anal, Chile, programmes 65.L-0428, 66.C-0446, Constraints for giant planetary companions to Proxima Cen us- 

267.C-5700, 68.C-0415, 69.C-0722, 70.C-0044, 71.C-0498, 072.C- ing HST Fine Guidance Sensor astrometry were reported by 
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precise RV data from the ESO Coude Echelle Spectrometer 
planet search were presented by Kiirster et al. (1999). The com- 
bination of both studies already excluded all companions with 
(minimum) masses higher than 0.8 Mjupitei for the period range 
1 to 600 days. 



2. Stellar properties of Proxima Cen 

Proxima Centauri (GJ 551, a Cen C, HIP 70890) is an M5Ve 
dwarf and with a distance of d=1.29 pc the closest star to 
the Sun. Hipparcos measured a parallax of 772.33 + 2.42 mas 
(Perryman et al. 1997). It is still under debate whether Proxima 
Cen is actually gravitationally bound to the a Cen AB binary 
(e.g. Wertheimer & Laughlin 2006). The star exhibits all char- 
acteristics of a magnetically active star, like coronal X-ray emis- 
sion (e.g. Hiinsch et al. 1998) and flare activity (e.g. Giidel et 
al. 2004). Benedict et al. (1998) reported on the photometric 
variability of Proxima Cen using 4 years of HST FGS photom- 
etry. They find a possible rotational period for Proxima Cen 
of x 84 days and an activity cycle of =s 1100 days. Recently 
Cincunegui et al. (2007) presented tentative evidence for a pos- 
sible activity cycle with a shorter period of ^ 442 days. 

A crucial parameter for the radial velocity method is the stel- 
lar primary mass. Using the V-band mass-luminosity relation- 
ship of Henry et al. (1999) we derive a mass of 0.108 M© for 
Proxima Cen, while the K-band mass-luminosity relationship of 
Delfosse et al. (2000) yields a slightly higher mass of 0.12 Mq. 



3. Observations and data reduction 

We observed Proxima Cen with the UVES spectrograph at the 
ESO VLT-UT2 as part of our ongoing Doppler search for low 
mass planets around a sample of 40 M dwarfs. We use a slit 
width of 0.3 arcsec and image slicer #3 to obtain a resolving 
power of R ^ 100, 000. The 37 orders on the two CCDs cover a 
wavelength range from 4950 to 7040 A. 

During 76 nights between March 2000 to March 2007 we 
collected a total of 229 individual spectra of Proxima Cen typi- 
cally grouped into three consecutive spectra per night. An expo- 
sure time of 780 seconds was selected for all spectra. The data 
have an average signal-to-noise ratio of53.1 +7.7 (per pixel) . 

To allow a simultaneous wavelength calibration, as well as 
the reconstruction of the instrumental profile, we employ the 
standard iodine (I2) cell technique for high precision RV mea- 
surements (e.g. Butler et al. 1996). The I2 spectrum has useful 
reference lines from x 5000 to 6500 A that are superimposed on 
the stellar spectrum. 

After bias subtraction, flat-fielding, subtraction of inter-order 
Echelle background, wavelength calibration, and barycentric 
correction using the JPL ephemeris DE200 (Standish 1990), we 
determine differential RVs with oni Austral code as described in 
Endl et al. (2000). 



4. RV results 

One result of the data modeling procedure is a formal internal 
uncertainty of the RV measurement. This internal uncertainty 
includes errors due to photon noise, CCD readout noise and "al- 
gorithmic" noise in the modeling process to remove all instru- 
mental effects (like spatial and temporal variations in the instru- 
mental profile). For the Proxima Cen data we obtained a mean 
internal error of 2.34 + 0.28 m s"' . 
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Fig. 1. 7 years of RV measurements for Proxima Cen using 
UVES+I2 cell at the ESO VLT. (The secular acceleration has 
been subtracted.) The data have a total rms of 3. 1 1 m s"' and an 
average uncertainty of 2.34 m s"' . 



As discussed in Kiirster et al. (2003) (from hereon called pa- 
per I) we then formed nightly means of the RV measurements 
from the typically three (but occasionally more) consecutive ex- 
posures. Removing this "dumpiness" of the data allows for are- 
liable estimation of the significance of possible signals using the 
bootstrap randomization method (see section 6). This reduced 
the total number of data points from 229 to 76. We then sub- 
tracted the expected secular RV acceleration of 0.45 ms"' yr ', 
due to the change in perspective to this nearby star, from the data 
(cf. paper I). Figure 1 displays our 76 differential RV measure- 
ments for Proxima Cen. The velocities are listed in Table [T] The 
data have a total rms of 3.1 1 m s"', slightly larger than the mea- 
surement uncertainties. The excess RV scatter is thus 2.05 m s"' . 

We noted that over the past 2 years the RV scatter of Proxima 
is greatly reduced. While the data taken before 2005 have an rms 
scatter of 3.04 m s"', the 16 data points we obtained since then 
have a scatter of only 1.36 ms"'. Since neither the instrument 
nor our data reduction algorithm was changed over this time, 
this change is likely instrinsic to the star 



4.1. Period search 

We use the classic Lomb-Scargle periodogram (Lomb 1976 ; 
Scargle 1982) to search for possible periodicities in the Proxima 
Cen RV data. The resulting power spectrum is shown in Fig. 
2. The highest peak is located at 363.6 days. There is also an 
intriguing peak at longer periods from 2000 to 3000 days, but 
such a long period is comparable to our observing time span. 
We find no power at the suspected rotation period of s; 84 d, 
and only moderate power at the possible activity cycle period of 
^ 1100 d. To estimate the significance of these signals we use 
the bootstrap randomization method (e.g. Kiirster et al. 1997). 
After 10 000 bootstrap randomizations of our data we find that 
the 364 d peak has a formal false-alarm-probability (FAP) of just 
0.05%. However, this peak is likely produced by the interplay of 
the strong 1-year peak in the window function (lower panel in 
Fig. 2) and the dense data cluster in the middle of 2001. This 
cluster is systematically lower than the rest of the points. If we 
reduce the weight of these points by averaging them we find a 
reduction in power of the 364-d peak. (The 1 -yr peak disappears 
completely if we remove the entire cluster from the time series.) 
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Fig. 2. Lomb-Scargle periodogram of the Proxima Cen RV data. 
The upper panel shows the power spectrum and the lower panel 
displays the window function. The 1-year peak in our window 
function is quite dominant. Using the bootstrap method to esti- 
mate the significance of the peak shows that no significant signal, 
except the peak close to 1 year, is present in our data. 



We will discuss the cause of the systematic blue-shift of these 
points in the next section. No other peak in the periodogram has 
a FAP less than 1%. 



5. Ha line index measurements 

As Proxima Cen is a known active flare star we suspected that 
the excess scatter might be caused by stellar activity and that we 
would find a correlation of the RV data with an activity indicator. 
In the same fashion as for Barnard's star (see paper I for details) 
we determined a line strength index for Hq, and, as a comparison, 
an index for a nearby Cal line whose strength should not depend 
on the activity level of the star. While for a relatively inactive 
M dwarf like Barnard's star the chromospheric emission is seen 
mostly as a small "filling in" of the H^ line core, in the case of 
Proxima Cen it is seen at all times as a strong emission feature 
above the continuum. We calculate the line index as the flux in- 
side a region of ±15 km s"' of the center of the line normalized 
to the flux of 2 adjacent spectral regions. As a check against sys- 
tematic errors we do the same for the Cal line at 6572.795 A. We 
determine the line index for each individual spectrum and then 
average the values for each night in the same manner as for the 
RV results. 

Proxima Cen reveals a large amount of variability in H„. The 
average H„ line index is 2.89 with a standard deviation of 0.90 
(31%), while the Cal indices is indeed quite constant. The mean 
Cal index is 0.57165 with a standard deviation of 0.0103 (1.8%). 
The amount of variability of the Hq. index is more than 4 times 
larger than what we have found in Barnard's star (7.6%). 

However, a period search in the line index data did not find a 
significant signal, likely due to the stochastic nature of the vari- 
ations (flaring). We find the highest power of 7. 1 3 at a period of 
25.7 days with an estimated FAP of 35%. 

Fig. [3] displays the Hq, line indices as a function of the RV 
results. We adopt as the error of the line index the rms of the 
individual measurements that were averaged. We find a weak 
correlation between the Hq, data and the RV results. The lin- 
ear correlation coefficient is -0.23 resulting in a 4.8% prob- 
ability that the null-hypothesis, that these two values are not 
correlated, is correct. A linear fit to the data yields a slope of 
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Fig. 3. Proxima Cen H,, line index measurements as a function 
of RV. The error bar of the line index is the rms of the individ- 
ual measurements (within one night) binned to form one data 
point, in the same manner as the RV data. There appears to be 
a weak correlation between the Hq, index and the RV data. The 
linear correlation coefficient is -0.23. A fit to the data is shown 
as a dashed line. The residual RV scatter around this slope is 
3.03 ms-'. 



-0.71 + 0.30 ms ' H„ - index '. The scatter around this fit is 
3.03 ms"' (as compared to the overall scatter of 3.11 ms"'). At 
least in a qualitative way this confirms the trend we have found 
in paper 1: spectra with a higher H„ index appear blue shifted. 
However, the correlation for Proxima Cen is much weaker and 
the slope is much less significant than in the case of Barnard's 
star. (See the detailed discussion in paper 1 why this might indi- 
cate a general convective redshift for M dwarfs.) It is interesting 
to note that despite the fact that Proxima Cen is a much more 
active star, its activity level does not produce a higher RV excess 
scatter ("jitter") than in Barnard's star 

Proxima Cen appears to have been slightly more active dur- 
ing the intensive campaign in summer 2001 . The average H^^ line 
index of this data cluster (16 points) is 3.33 + 0.29 (error of the 
mean) while the mean of the rest of the data is 2.77 + 0.1 1 (60 
points). This small (1.4cr) difference might be the explanation 
why this group of measurements is systematically blue-shifted 
from the rest. 



6. Limits for planetary companions 

With the absence of a clear RV signal of planets in the Proxima 
Cen data we pose the question which type of planets we could 
have already detected and hence we can exclude. We compute 
mass upper limits with a method of injecting and recovering sig- 
nals into the original RV data set (Endl et al. 2001). We start 
with the assumption that the residual scatter in our data is the 
best representation of the errors and use the bootstrap random- 
ization method to determine the statistical significance of the re- 
covered signal. We intentionally do not correct the RV data with 
the small activity related slope we determined in the previous 
section, in order to determine our detection efficiency also in the 
presence of this source of correlated noise. The mass limit for 
a given period is set at the lowest signal amplitude where the 
signal is recovered at all trial phase angles with a significance 
of > 99%. The test signals were generated at 16 different phase 
angles, each shifted by 7r/8. We restrict the simulations to circu- 
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lar orbits. Limits derived for eccentric orbits are typically higher 
(see e.g. Wittenmyer et al. 2006 for a comparison between e = 
and e < 0.6 limits), but in this case we are particularly interested 
in planets that reside inside the habitable zone (HZ) all the time. 
An additional constraint is that for a successful recovery of the 
signal the peak in the power spectrum has to be at (or near) the 
correct input period. This is especially important for Proxima 
where the 1 -yr period at the maximum of the window function 
appears frequently as the dominant signal. In a conservative ap- 
proach we adopt the higher value of 0. 12 M© for the stellar mass 
(see section 2). 

FiglH shows the mass upper limits for planets based on this 
method. The solid line represents the mass range where all test 
signals were recovered at (or near) the correct period with a sig- 
nificance level of > 99%. The shaded area delimits the classic 
liquid surface water HZ after Kasting et al. (1993). For such a 
late type star the HZ is very close to the host star, for Proxima 
Cen it ranges from 0.022 to about 0.054 AU, corresponding to 
orbital periods from 3.6 to 13.8 days. 

We can not determine reliable upper limits in the period 
range close to 1 year. The fact that a strong signal at this pe- 
riod is present in the original data leads to the effect that even 
a zero amplitude input signal in phase with the original signal 
will always lead to a successful detection. We therefore exclude 
the region of 300 < P > 400 days from our simulations. (This 
excluded part of the mass-period diagram is shown as "1-year 
window" in FigH]) 

For periods outside the gap at 1 year we can exclude all plan- 
ets with msini > 16 M®, i.e. planets with minimum masses 
greater than Neptune's mass, out to 1 AU. For periods less than 
100 days (or a < 0.21 AU) we could have detected all "Super- 
Earths" with m sin ; > 8.5 M®). For the HZ of Proxima Cen we 
can rule out the presence of all planets with m sin /> 2 - 3 M® 
in circular orbits. In order not to confuse these upper limits with 
our actual sensitivity to low mass planets we also show the mass 
range where 75, 50 and 25% of the test signals were recovered 
with a statistical significance of > 99% (dashed, dash-dotted 
and long-dashed lines in Fig|4]i. This demonstrates that we al- 
ready had a 50% chance to detect planets with masses as low as 
m sin i - 1 .5 M® inside the HZ of Proxima Cen and slightly in- 
side the inner edge of the HZ even down to msini a; 1.2 M^. 
Over the entire range of periods (with the exception of the 
1-yr window) we are 100% (75, 50 and 25%) complete for 
RV signals with semi-amplitudes of K - AAA + 0.33 ms"' 
(3.47 + 0.32, 2.56 + 0.44 and 2.01 + 0.60 m s"', respectively). 

7. Discussion 

We began our VLT/UVES Doppler survey to search for terres- 
trial planets in the HZ of M dwarfs already in 2000. We pre- 
sented first estimates of the survey sensitivity for low mass plan- 
ets in paper I for Barnard's star and in Endl et al. (2003) for 
Proxima Cen. With the results from this paper we demonstrate 
that the discovery of m sin / ^ 1 M® is within our grasp. Since 
sensitivity is a function of RV precision, number of measure- 
ments and sampling, adding more points to the existing data 
string in a pseudo-random fashion, will allow us to improve the 
detection sensitivity over time. 

Guedes et al. (2008) discuss an RV search for terrestrial plan- 
ets in the HZ of a Cen B. They conclude that the detection of 
a 1 .7 Me planet will require continuous monitoring of the star 
over at least 3 years and the accumulation of nearly 100 000 in- 
dependent measurements with a precision of s; 3ms"'. We have 
shown here that by switching targets to Proxima Cen the same 



feat can be achieved with less than 100 data points. The main 
difference being the lower mass of the star as well as the much 
shorter orbital periods of HZ planets. 

When interpreting the mass limits, the reader should bear in 
mind that we only considered circular orbits. Limits for plan- 
ets on eccentric orbits are typically slightly higher (see e.g. 
Wittenmyer et al. 2006). Planets with masses above our mass 
threshold for circular orbits can still exist around Proxima Cen 
on eccentric orbits. We also considered only the case of a single 
planet. The RV signals of a multi-planet system with several low- 
mass bodies, is likely to be more difficult to detect by a pure pe- 
riodogram analysis (depending on their period spacing and mass 
ratios) and requires a significantly larger data sets. Simulations 
to determine the mass limits for multiple planets is beyond the 
scope of this paper 

The question whether terrestrial planets inside the HZ of M 
dwarfs are suitable for life has been discussed extensively in the 
literature (see e.g. Tarter et al. 2007 and Scalo et al. 2007 for 
a review of this subject). At first glance, M dwarf HZ planets 
did not appear to be attractive objects for astrobiology, because 
their small semi-major axes would lead to tidal locking into syn- 
chronous rotation and presumably to atmospheric collapse at the 
night side of the planet. However, more detailed climate stud- 
ies have shown that energy transfer between the day and night 
side of such a planet can prevent atmospheric collapse (Joshi et 
al. 1997). Furthermore, the long main-sequence life times of M 
dwarfs would presumably lead to stable environments on such 
planets on very long time scales. Recently, Lissauer (2007) and 
Raymond et al. (2007) discuss the potential of habitable M dwarf 
HZ planets in the context of their formation and in particular 
water delivery and find that these planets might be devoid of 
volatiles and thus less suitable for life. 

One of the drawbacks of the RV technique is the sin / am- 
biguity, we measure only minimum mass values while the true 
mass remains unknown. But mass is, of course, a crucial param- 
eter for the characterization of a planet. Low mass planets inside 
the HZ of M dwarfs have a higher probability to transit their host 
star, due to their smaller semimajor axes. Moreover, the transit 
depth for a given planet will be deeper for M dwarfs, than for 
earlier spectral types (see e.g. Gillon et al. 2007). Finding more 
short periodic transiting Super-Earths around M dwarfs will thus 
potentially increase the sample of planets for which we can dis- 
tinguish their internal composition by putting them on a mass- 
radius diagram (e.g. Selsis et al. 2007 ; Valenica et al. 2007). 

We have shown that finding terrestrial planets in short peri- 
ods around M dwarfs is possible even without an RV precision 
of 1 m s"' or better. M dwarfs with masses less than 0.2 M© are 
numerous in the solar neighborhood. However, the majority of 
them are too faint in the V-band, to be accessed by traditional 
Doppler surveys, working in the optical spectral range. A near 
infrared high resolution spectrograph (e.g. Guenther et al. 2006) 
that can obtain an RV precision of 2 - 3 m s ' would thus be the 
ideal tool to carry out a search for nearby m sin / < 1 M® planets. 

8. Conclusions 

1 . We present 7 years of high precision RV data for our clos- 
est neighbor in space, the M5V star Proxima Cen, obtained 
with UVES + I2 cell at the ESO VLT/UT2. We detect no 
significant periodicities (except close to the 1-yr peak in the 
window function) that can be attributed to orbiting compan- 
ions. 

2. Using the same set of spectra we measure an H„ line index to 
estimate the magnetic activity level of Proxima Cen. These 
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Fig. 4. Mass upper limits for planets in circular orbits around Proxima Cen based on our numerical simulations. All test signals 
with amplitudes corresponding to masses on and above the solid line were recovered with a > 99% significance. The solid, dashed, 
dashed-dotted and long-dashed lines show the mass range where we successfully recovered 100, 75, 50 and 25% of the test signals. 
The shaded area labelled "HZ" displays the approximate location of the classic liquid water habitable zone (after Kasting et al. 1993). 
The region labelled "1-year window" shows the period range that we excluded from our simulations. 



line indices show a large amount of scatter due to flaring 
activity and are only weakly correlated with the RV results. 
3. Based on numerical simulations we demonstrate that we 
could have already detected all planets with msini = 2 - 
3 Me on circular orbits inside the classic habitable zone of 
Proxima Cen (assuming a stellar mass of 0.12 M®). 
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Table 1. Differential radial velocities of Proxima Cen from VLT/UVES 



JD [2 400 000+] 


dRV[ms-'] 


cr [ms '] 


51634.7343 


7.86 


2.16 


51654.6206 


0.31 


2.38 


51681.7798 


-0.28 


2.66 


51707.6289 


-0.34 


2.17 


51737.5345 


-0.26 


3.23 


51740.5824 


3.34 


2.22 


51761.5579 


-1.96 


2.11 


51818.4977 


-1.79 


2.52 


51919.8414 


1.55 


2.18 


51946.8066 


2.47 


2.05 


52118.5200 


-7.35 


2.59 


52119.5002 


-4.72 


2.57 


52120.4995 


-4.52 


2.19 


52121.5046 


-1.50 


2.21 


52122.5041 


-2.27 


2.23 


52123.5158 


-0.99 


2.24 


52124.5091 


-1.62 


2.26 


52125.5043 


-1.92 


2.20 


52126.5086 


-1.09 


2.13 


52127.5023 


-3.77 


2.16 


52128.5090 


-2.84 


2.24 


52129.5172 


-2.64 


2.19 


52130.5194 


-5.14 


2.05 


52131.5036 


-6.23 


2.11 


52134.5201 


-1.61 


2.21 


52135.5366 


-1.00 


2.31 


52308.8586 


-0.16 


2.00 


52320.8412 


2.17 


2.38 


52336.8520 


3.05 


2.19 


52380.6817 


-1.15 


2.30 


52394.7349 


-3.12 


2.15 


52415.6011 


-0.46 


2.85 


52424.6250 


1.45 


2.23 


52425.6540 


-2.85 


2.73 


52437.5534 


-1.41 


2.06 


52482.6464 


-2.09 


2.91 


52490.5975 


-1.93 


2.77 


52530.5056 


-5.08 


2.54 


52532.5082 


-0.84 


2.21 


52657.8318 


3.25 


2.25 


52659.7972 


1.70 


2.32 


52684.8259 


3.01 


2.56 


52709.6935 


4.08 


2.14 


52743.8400 


-0.97 


2.10 


52761.5228 


0.32 


2.11 


52767.9010 


-5.19 


2.93 


52788.5002 


0.49 


2.01 


53071.8829 


-1.35 


2.18 


53076.8916 


-5.27 


2.30 


53088.8399 


1.37 


2.07 


53090.8887 


-0.67 


2.40 


53095.8975 


4.80 


2.37 


53096.8803 


3.81 


3.49 


53098.7983 


0.76 


2.18 


53112.6938 


-5.49 


2.48 


53118.5595 


3.86 


2.46 


53125.6977 


-1.07 


2.40 


53128.6031 


3.79 


2.22 


53133.5396 


-1.94 


2.67 


53207.4883 


-2.18 


2.50 
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Table 1. continued. 



JD [2 400 000+] 


dRV[ms-'] 


cr [ms '] 


53394.8511 


0.26 


2.17 


53413.8344 


3.13 


2.32 


53430.8431 


1.76 


2.11 


53443.8731 


4.69 


2.23 


53453.7691 


4.56 


2.12 


53465.8306 


2.74 


2.64 


53487.6901 


2.89 


2.18 


53495.7335 


0.61 


2.05 


53520.5247 


2.50 


2.92 


53748.8611 


3.90 


2.45 


53783.7137 


3.88 


2.50 


53813.6861 


1.40 


2.18 


53823.6115 


2.51 


2.33 


54170.8614 


4.20 


2.15 


54188.8877 


1.40 


2.29 


54189.7161 


3.12 


2.32 



